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Abstract 

We  consider  the  problem  of  approximating  normal  and 
feature  sizes  of  a  surface  from  point  cloud  data  that  may 
be  noisy.  These  problems  are  central  to  many  applica¬ 
tions  dealing  with  point  cloud  data.  In  the  noise-free 
case,  the  normals  and  feature  sizes  can  be  approximated 
by  the  centers  of  a  set  of  unique  large  Delaunay  balls 
called  polar  balls.  In  presence  of  noise,  polar  balls  do 
not  necessarily  remain  large  and  hence  their  centers  may 
not  be  good  for  normal  and  feature  size  approximations. 
Earlier  works  suggest  that  some  large  Delaunay  balls  can 
play  the  role  of  polar  balls.  However,  these  results  were 
short  in  explaining  how  the  big  Delaunay  balls  should  be 
chosen  for  reliable  approximations  and  how  the  approx¬ 
imation  error  depends  on  various  factors.  We  provide 
new  analyses  that  fill  these  gaps.  In  particular,  they  lead 
to  new  algorithms  for  practical  and  reliable  normal  and 
feature  approximations. 

1.  Introduction 

Recently,  a  number  of  algorithms  have  been  designed 
for  processing  point  cloud  data.  Often  these  algorithms, 
as  a  basic  step,  estimate  the  normals  and  features  of 
the  sampled  surface  from  the  given  point  cloud.  For 
example,  some  algorithms  [1,  8,  11]  need  a  normal 
estimation  step  for  surface  reconstruction,  and  others 
estimate  the  scale  of  local  geometry  also  called  the  local 
feature  size  to  handle  non-uniform  samples  [9,  14].  In 
the  noise-free  case  the  problem  of  normal  and  feature 
size  approximations  have  been  well  studied  [2,  4,  6].  In 
the  case  of  noise,  optimization  based  techniques  [1,  13] 
are  known  for  normal  approximations  though  they  do 
not  have  theoretical  guarantees.  It  is  known  that  results 
from  the  noise-free  case  can  be  extended  by  using  big 
Delaunay  balls  that  can  help  in  estimating  normals  [8] 
with  theoretical  guarantees.  However,  it  is  not  known 
how  the  error  of  approximation  depends  on  different 
noise  components,  and  more  importantly,  how  the  big 
Delaunay  balls  should  be  chosen  for  reliable  approx¬ 
imations.  The  problem  for  feature  approximations  in 
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presence  of  noise  is  much  less  understood.  No  reliable 
and  practical  algorithm  is  known  for  it.  In  this  paper  we 
address  these  open  issues. 

Motivation  and  results.  Amenta  and  Bern  [2]  intro¬ 
duced  the  concept  of  poles.  These  are  the  furthest 
Voronoi  vertices  from  the  respective  sites  on  the  two 
sides  of  the  sampled  surface.  In  terms  of  the  Delaunay 
triangulation,  poles  are  the  centers  of  the  largest  Delau¬ 
nay  balls  incident  to  the  sample  points  on  both  sides  of 
the  sampled  surface.  These  balls  are  also  called  the  po¬ 
lar  balls.  Amenta  and  Bern  showed  that,  in  the  noise-free 
case,  the  normals  can  be  estimated  by  the  poles.  Further, 
Amenta,  Choi,  Kolluri  [4]  and  Boissonnat,  Cazals  [6] 
proved  that  the  poles  also  approximate  the  medial  axis 
and  hence  local  feature  sizes  can  be  estimated  by  dis¬ 
tances  to  the  poles. 

In  the  presence  of  noise,  the  above  results  do  not  hold 
since  some  of  the  polar  balls  can  be  arbitrarily  small 
with  their  centers  being  arbitrarily  close  to  the  surface. 
See  the  top-middle  picture  in  Figure  1.  Nevertheless, 
Dey  and  Goswami  [8]  observed  that,  under  a  reasonable 
noise  model,  many  Delaunay  balls  remain  big  and  their 
centers  can  help  in  approximating  the  normals.  The  error 
in  the  normal  approximation  by  big  Delaunay  balls  ob¬ 
viously  depends  on  the  sampling  density  ( e )  and  also  on 
the  size  of  the  chosen  Delaunay  balls.  A  detailed  anal¬ 
ysis  on  these  dependencies  is  missing  in  earlier  works. 
First,  our  analysis  provides  an  error  bound  that  unifies 
earlier  results.  Second,  it  tells  us  how  to  choose  big  De¬ 
launay  balls  in  practice  for  reliable  normal  and,  in  partic¬ 
ular,  feature  size  approximations  for  noisy  point  clouds. 

In  the  noise-free  case  the  choice  of  the  Delaunay  balls 
is  not  an  issue  in  normal  and  feature  size  approxima¬ 
tions  as  they  are  approximated  from  polar  balls  which 
are  almost  as  big  as  medial  balls.  However,  in  the  case 
of  noise,  the  choice  of  Delaunay  balls  is  an  issue  as  all 
polar  balls  are  not  big.  To  remain  scale  independent  one 
can  choose  Delaunay  balls  whose  radii  are  larger  than 
a  threshold  determined  by  some  nearest  neighbor  dis¬ 
tances  of  the  sample  points  incident  on  the  Delaunay 
balls.  In  order  to  gauge  the  viable  values  of  the  thresh¬ 
old,  it  is  important  to  know  how  the  normal  and  feature 
size  approximation  errors  depend  on  the  radii  of  the  De¬ 
launay  balls.  Our  new  analysis  provides  this  relation. 
We  show  that  normals  can  be  estimated  from  Delaunay 
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Figure  1:  Top  row:  Left:  noise-free  case,  poles  are  approximating  the  medial  axis  and  normals  well.  Middle:  A  small 
noise  disturbs  the  poles  significantly  resulting  in  poor  normal  and  medial  axis  approximation  with  all  poles.  Right: 
only  a  subset  of  big  Delaunay  balls  are  chosen,  normals  though  not  medial  axis  are  well  approximated.  Bottom  row: 
Left:  Delaunay  balls  of  bigger  size  are  chosen  to  exclude  unwanted  poles,  some  significant  parts  of  the  medial  axis  are 
not  approximated.  Right:  Centers  of  polar  balls  chosen  with  our  algorithm  approximate  the  medial  axis  everywhere. 
Approximated  feature  sizes  are  indicated  in  the  highlighted  boxes. 


balls  that  are  not  necessarily  as  big  as  local  feature  sizes 
(/(•)).  In  fact,  Delaunay  balls  with  radii  as  small  as 
£2  /(•)  are  also  good  for  normal  estimations.  See  top 
row  in  Figure  1  for  an  illustration. 

The  case  for  feature  estimations  in  presence  of  noise 
is  far  more  difficult.  This  is  because,  unlike  normal  ap¬ 
proximations,  not  all  centers  of  Delaunay  balls  chosen 
with  a  reasonable  threshold  approximate  the  medial  axis. 
Choosing  the  right  ones  is  hard.  If  the  threshold  is  rel¬ 
atively  small,  a  number  of  centers  remain  which  do  not 
approximate  the  medial  axis.  See  top-right  picture  in 
Figure  1 .  On  the  other  hand  if  the  threshold  is  large,  the 
medial  axis  for  some  parts  of  the  models  may  not  be  ap¬ 
proximated  at  all;  see  bottom-left  picture  in  Figure  1 .  As 
a  result  no  threshold  may  exist  for  which  large  Delaunay 
balls’  centers  approximate  the  medial  axis,  the  DOG  data 
in  Figure  1  and  the  HORSE  data  in  Figure  8  are  two  such 


examples  in  two  and  three  dimensions  respectively. 

We  propose  a  different  algorithm  to  choose  the 
Delaunay  balls  for  approximating  the  medial  axis.  We 
consider  ^-nearest  neighbors  for  some  k  and  take  the 
largest  polar  ball’s  center  among  these  neighbors  to 
approximate  the  medial  axis.  Our  analysis  leads  to  this 
algorithm  which  frees  us  from  the  burden  of  choosing 
a  size  threshold.  Our  experiments  suggest  that  k  can 
be  chosen  fairly  easily,  generally  in  the  range  of  5  to 
10.  The  most  important  thing  is  that  a  k  can  be  found 
for  which  the  medial  axis  is  well  approximated  where 
no  such  size  threshold  may  exist.  The  bottom  row  of 
Figure  1  illustrates  this  point. 

Previous  results.  Amenta,  Bern  and  Eppstein  [3] 
introduced  the  e-sampling  for  noise-free  case.  This 
requires  each  point  on  the  surface  to  have  a  sample 
point  within  a  distance  of  e  times  the  local  feature  size. 
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When  noise  is  allowed,  the  sample  points  need  not 
lie  exactly  on  the  surface  and  may  scatter  around  it. 
Therefore,  the  sampling  model  needs  to  specify  both  a 
tangential  scatter ,  i.e.,  the  sparseness  of  the  sampling 
along  tangent  directions  of  the  surface  and  also  a  normal 
scatter ,  i.e.,  the  sparseness  of  sampling  along  the  normal 
directions.  Dey  and  Goswami  [8]  introduced  a  noise 
model  that  uses  the  same  sampling  parameter  e  for  both 
scatters.  Kolluri  [11]  and  later  Dey  and  Sun  [9]  modified 
the  normal  scatter  to  have  £ 2  dependence.  The  errors 
of  normal  and  feature  approximations  depend  on  both 
tangential  and  normal  scatters.  Therefore,  we  introduce 
two  independent  parameters  e  and  5  for  these  two 
scatters  to  reveal  the  dependence  of  the  approximation 
errors  on  these  two  parameters  separately. 

Normal  approximation:  Dey  and  Goswami  [8]  and 
Mederos  et  al.  [12]  showed  that  when  both  tangential 
and  normal  scatters  are  0(e)  times  the  local  feature  size, 
the  normals  can  be  approximated  with  an  0(\/^)  error  if 
the  chosen  Delaunay  balls  have  radius  almost  as  big  as 
the  local  feature  size.  Dey  and  Sun  [9]  showed  that  the 
error  is  0(e)  if  the  normal  scatter  is  only  0(£2)  times 
the  local  feature  size.  None  of  these  results  specify  how 
the  error  depends  on  the  radii  of  the  chosen  Delaunay 
balls. 

In  this  paper  we  provide  a  simple  elegant  analysis 
which  shows  that  the  error  is  2(A  +  1  )0(£  +  s/d) 
where  A  is  the  radius  of  the  Delaunay  ball.  Previous 
results  under  different  noise  models  can  be  derived 
from  this  unified  result.  One  implication  of  this  result 
is  that  Delaunay  balls  as  small  as  0(^2  +  5^)f^)  can 
help  in  estimating  the  normals.  This  relaxes  the  burden 
on  setting  the  parameter  for  the  normal  estimation 
algorithm. 

Feature  approximation:  Amenta,  Bern  and  Epp- 
stein  [3]  defined  the  local  feature  size  of  a  point  x  on 
the  surface  as  the  distance  of  x  to  the  medial  axis.  Ob¬ 
viously,  the  local  feature  size  can  be  estimated  if  the 
medial  axis  can  be  approximated.  An  algorithm  for  ap¬ 
proximating  the  medial  axis  from  noisy  point  clouds  ex¬ 
ists  [7].  This  algorithm  approximates  the  medial  axis 
with  Voronoi  faces  under  a  stringent  uniform  sampling 
condition.  Selecting  Voronoi  faces  to  approximate  the 
medial  axis  is  not  a  simple  task  in  practice  even  for 
noise-free  case  [5,  10]  and  it  is  not  clear  how  this  al¬ 
gorithm  works  in  practice  when  noise  is  present.  More¬ 
over,  for  estimating  the  local  feature  size  a  continuous 
approximation  with  Voronoi  faces  is  an  overkill.  A  dis¬ 
crete  approximation  of  the  medial  axis  with  a  set  of 
Voronoi  vertices  serves  the  purpose  equally  well.  For 
the  noise-free  case,  such  an  approximation  was  proposed 
by  Amenta  et  al.  [4]  and  Boissonnat  and  Cazals  [6].  Re¬ 


cently,  Mederos  et  al.  [12]  derived  some  results  for  noisy 
point  clouds  that  have  some  connections  to  the  local  fea¬ 
ture  size  approximations  though  the  approximation  fac¬ 
tor  depends  on  a  surface  related  constant  which  can  be 
potentially  huge. 

Our  analysis  is  free  of  any  surface  dependent  constant 
and  it  relates  the  approximation  error  to  the  tangential 
and  normal  scatters  separately.  Most  importantly,  the 
analysis  justifies  our  choice  of  polar  balls  based  on  near¬ 
est  neighbors  to  approximate  the  medial  axis.  Figure  1 
and  8  show  that  this  choice  is  far  more  superior  than  the 
big  Delaunay  ball  strategy.  Experiments  with  our  imple¬ 
mentation  [15]  of  the  algorithm  confirm  this  claim  for 
other  models. 

2.  Preliminaries 

2.1  Definitions 

For  a  set  Y  C  R3  and  a  point  x  G  R3,  let  d(x.  Y)  denote 
the  Euclidean  distance  of  x  from  Y ;  that  is, 

d(x,Y)  =  inf  {||j/- *||). 
y£Y 

The  set  Bc,r  =  {y  \  y  G  M3,  \  \y  —  c\\  <  r}  is  a  ball  with 
radius  r  and  center  c. 

Voronoi  and  Delaunay  diagram.  For  a  finite  point  set 
P  C  M3,  we  will  denote  the  Voronoi  diagram  and  its 
dual  Delaunay  triangulation  of  P  by  Vor  P  and  Del  P 
respectively.  The  Voronoi  cell  for  a  point  p  is  denoted  as 

Vp. 

Sampled  surface.  Let  S  C  l3  be  a  compact  smooth 
surface  without  boundary  from  which  the  input  sample 
is  derived  possibly  with  noise.  Also,  assume  that  E  is 
connected.  The  bounded  and  unbounded  components  of 
M3  \  E  are  denoted  Dj  and  Do  respectively.  The  nor¬ 
mal  at  any  point  x  E  E  is  denoted  which  is  directed 
locally  inward,  i.e.,  toward  Dj. 

The  medial  axis  M  of  E  is  the  locus  of  the  centers  of 
the  maximal  balls  whose  interiors  are  empty  of  points 
in  E.  These  balls  meet  E  only  tangentially.  We  call 
each  such  ball  Bm^r  a  medial  ball  where  r  =  d(m ,  E). 
Barring  some  pathological  cases,  we  can  assume  Mr I  E 
is  empty  if  E  is  smooth.  The  subsets  of  M  in  Dj  and 
Do  are  called  inner  and  outer  medial  axis  respectively. 
For  each  point  x  G  E,  there  are  two  medial  balls ,  one 
centering  a  point  in  the  inner  medial  axis  and  the  other 
in  the  outer  medial  axis.  The  local  feature  size  at  a  point 
x  e  E  is  defined  as  f(x)  =  d(x ,  M).  The  function  /() 
satisfies  the  following  Lipschitz  property  [2]. 

Lipschitz property.  For  any  two  points  x,  y  G  E,  f(x)  < 

f  ly)  +  Ik  -  y II- 
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2.2  Sampling 

A  finite  set  of  points  P  C  Sis  called  an  e-  sample  of  £  if 
d(Xj  P)  <  ed(Xj  M)  for  each  x  £  £.  To  accommodate 
the  tangential  and  normal  scatters  of  points  around  £  in 
the  noisy  case,  we  put  two  conditions  on  the  sampling. 
The  first  condition  says  that  the  projection  of  the  point 
set  P  on  the  surface  makes  a  dense  sample  and  the 
second  one  says  that  P  is  close  to  the  surface.  We 
also  use  a  third  condition  to  make  the  sampling  locally 
uniform.  To  make  the  sampling  definition  general,  we 
use  a  separate  parameter  for  each  sampling  condition. 
For  any  point  x  £  M3  \  M  let  x  denote  its  closest  point 
on  £.  Clearly,  the  segment  xx  is  parallel  to  the  normal 
n5. 

We  say  P  C  M3  is  a  (e,  5,  k) -sample  of  £  if  the  follow¬ 
ing  conditions  hold. 

(i)  P  =  {p}pep  is  an  6-sample  of  £, 

(ii)  \\p-p\\  <  Sf(p), 

(iii)  | |p  —  q\ |  >  ef  (p)  for  any  two  points  p,  q  in  P  where 
q  is  the  ftth  nearest  sample  point  to  p. 

Figure  2  illustrates  why  we  put  the  third  condition.  In 
the  figure  the  same  point  sample  satisfies  the  first  two 
conditions  for  two  different  curves;  C  and  also  C  U  Cf . 
Our  analyses  for  normal  and  medial  axis  approximations 
apply  to  both  of  these  curves;  albeit  with  different  scales 
of  local  feature  sizes.  Therefore,  the  analyses  do  not 
need  the  third  condition  in  the  sampling.  However,  our 
approximation  algorithms  determine  a  particular  scale 
by  looking  at  the  nearest  neighbor  distances.  This  im¬ 
plies  that  the  sampling  cannot  allow  the  ambiguity  which 
is  forced  by  assuming  a  local  uniformity  constraint  in  the 
third  one. 


Figure  2:  A  point  sample  satisfying  sampling  conditions 
(i)  and  (ii)  for  a  single  component  curve  C  (left)  and  also 
the  curve  CuC"  (right). 

In  the  analysis  we  concentrate  only  in  the  bounded 
component  Qi  together  with  the  inward  normals  and  in¬ 
ner  medial  axis.  It  should  be  clear  that  the  results  also 
hold  for  unbounded  component,  outward  normals  and 
outer  medial  axis.  For  a  point  x  £  £,  let  mx  denote  the 
center  of  the  inner  medial  ball  meeting  £  at  x  and  px  its 
radius. 

It  follows  almost  immediately  from  our  sampling  con¬ 
ditions  that  each  point  of  E  and  a  point  not  far  away  from 


E  has  a  sample  point  nearby.  Lemma  1  and  Corollary  2 
formalize  this  idea. 

Lemma  1.  Any  point  x  £  E  has  a  sample  point  within 
Sif(x)  distance  where  s\  —  (6  +  e  +  &). 

Proof.  From  the  sampling  condition  (i),  we  must  have  a 
sample  point p  so  that  || x—p\ |  <  ef(x).  Also,  \\p— p\\  < 
Sf(p)  <  <5(1  +  e)f(x).  Thus, 

Ik-Pll  <  \\x  ~  P\\  +  \\P  ~  p\\ 

<  ef(x)+S(l  +  e)f(x) 

<  (S  +  e  +  6e)f(x). 

□ 

Since  f(x)  <  px  for  any  point  x  €  S,  the  following 
corollary  is  immediate. 

Corollary  2.  Any  point  y  Gl3  with  \\y  —  y\\  =  6  py  has 
a  sample  point  within  62 py  distance  where  62  =  (25  + 
e  +  &). 

3.  Empty  balls 

A  ball  is  empty  if  its  interior  is  empty  of  points  from  P. 

A  main  ingredient  in  our  analysis  will  be  the  existence  of 
large  empty  balls.  They  in  turn  lead  to  the  existence  of 
large  Delaunay  balls  that  circumscribe  Delaunay  tetrahe- 
dra  in  Del  P.  The  centers  of  such  Delaunay  balls  which 
are  also  Voronoi  vertices  in  Vor  P  play  crucial  roles  in 
the  algorithms  for  normal  and  feature  estimations.  In 
this  section,  we  present  two  lemmas  that  assure  the  exis¬ 
tence  of  large  empty  balls  with  certain  conditions. 

Lemma  3  below  assures  that  for  each  point  x  G  E 
there  is  a  large  empty  ball  of  radius  almost  as  large  as 
(i)  f(x)  and  (ii)  px.  Notice  the  differences  between  the 
distances  of  these  balls  from  x.  Also,  see  Figure  3. 

Lemma  3.  A  ball  is  empty  of  sample  points  from 
P  if  either 

(i)  fri  =  x,  || m  —  x\\  =  f(x)  and  r  =  (1  —  3 S)f(x), 
or 

(ii)  m  =  mx  and  r  =  (1  —  S)px. 

Proof  Let  p  be  any  point  in  P  (the  left  picture  of  Fig¬ 
ure  3).  For  (i)  we  have 

f(P)  <  f(x)  +  \\x-p\\  <  f(x) +  \\x -m\\  +  \\m-p\\ 

<  2f(x) +  \\m-p\\. 

Then, 

\\m-p\\  >  \\m-p\\-\\p-p\\>\\m-p\\-5f(p) 

>  \\m-p\\-S(2f(x)  +  \\m-p\\) 

=  (l-5)\\m-p\\-2Sf(x) 

>  (1-3  5)f(x) 
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Figure  3:  Illustration  for  Lemma  3.  The  dotted  big  balls 
are  not  empty  of  sample  points  but  their  slightly  shrunk 
copies  (shown  with  solid  boundaries)  are. 


Figure  4:  Illustration  for  Lemma  4. 


as  ||ra— pH  >  ||ra  —  x\\  =  f(x).  Therefore,  p  cannot  be 
in  the  interior  of  B.  Now  consider  (ii).  We  get 


l|™*-p||  > 
> 

> 

> 


\\mx  -pH  -  lb -pH 
Ib^  -P\\-f>f(p) 

II mx  -p\\  -5\\mx  -p\\ 
(1  -5)\\mx  -p|| 

(1  -  S)Px 


as  llm*  -p||  >  ||mx  -  a:||  =  px. 


□ 


Next,  we  show  that,  for  each  point  x  of  S,  there  is 
a  nearby  large  ball  which  is  not  only  empty  but  also 
its  boundary  passes  through  a  sample  point  close  to  x. 
Eventually  these  balls  will  be  deformed  to  Delaunay 
balls  for  medial  axis  approximations. 

Lemma  4.  For  each  point  x  £  E  there  is  an  empty  ball 
Bc^r  with  c  £  Ft i  that  enjoys  the  following  properties: 

(i)  r  is  at  least  (1  —  2v/£^)piC,  mx  is  in  Bc^r  and  ||c  — 
mx\\  <  2 ydFjpx  where  e 2  defined  in  Corollary  2  is 
0(s  +  5), 

( ii )  its  boundary  contains  a  sample  point  p  within  a  dis¬ 
tance  63 px  from  x  where  63  =  2e2  ¥  +  5  and  e,  5 
are  sufficiently  small. 


Proof  Consider  the  empty  ball  B  =  Bmx whose 
boundary  passes  through  a  point  y  where  y  =  x  and 
\\y  —  x\\  =  5px  and  R  =  (1  —  5)px.  Such  a  ball  exists 
by  Lemma  3  (ii)  (see  Figure  4). 

Shrinking:  For  a  ft  <  1,  let  B 13  =  Bmx^R  be  a 
further  shrunk  copy  of  B.  The  ball  B  and  hence  B@  are 
empty. 

Rigid  motion:  Translate  B$  rigidly  by  moving  the 
center  along  the  direction  vrCxx  until  its  boundary  hits  a 


sample  point  p  £  P.  Let  this  new  ball  be  denoted  Bc,r, 
refer  to  Figure  4. 

Obviously  r  =  f3R.  Let  d  =  \\c— mx\\.  First  we  claim 
(1  —  f3)R  <  d  <  (1  —  fi)R  +  62 Px- 


The  first  half  of  the  inequality  holds  since  B  is  empty 
of  samples  and  hence  B@  has  to  move  out  of  it  to  hit  a 
sample  point.  The  second  half  of  the  inequality  holds 
since  from  Corollary  2,  a  ball  centered  at  y  with  radius 
£2 px  cannot  be  empty  of  samples. 

Next  we  obtain  an  upper  bound  on  \\y  —  p\\.  Refer  to 
Figure  4.  Since  \\p'  —  q\\2  =  ||c  —  p'||2  —  ||c  —  q\\2  = 
\\mx  ~  p'\ |2  —  ||ra^  —  q\\2,  we  have 

\\c-p,\\2-\\c-q\\2  =  R2-(\\c-q\\  +  d)2 


which  gives 


c-bl 


R2_r2_cp 

2d 


Hence 


Ib-Pll2 


< 


< 


lb-p'll2 

lb'-bl2  +  lk-p|l2 

tf-(d+\\c-q\\f  +  (R-(d+\\c-q\\)f 


2 R2  -Rd-  ^(R2  -r2) 

£2(j-  +P)  p2 
{l-8){l-fi)F£2  ’ 


(3.1) 


Since  we  want  both  1 1  c  —  mx  \  |  and  1 1  x  —  p \  |  to  be  small, 
we  take  fi  —  1  —  2 .  Then  r  =  (3R  >  (1  —  2^/ef)px. 

Also  ||c  -  ra^H  <  ((1  -  5)^  +  £2)px  <  2y^2 Px, 
which  is  small  compared  to  r  given  e  and  5  are  suffi¬ 
ciently  small.  So  mx  stays  inside  Bc^r.  In  addition,  from 
inequality  3.1  we  have 


\\y-p\\  < 


(  2~v^ 

\1  —  <5  + 


&2  —  ^2  Px  • 


The  bound  on  \\p  —  x\\  follows  as  \\x  —  y\\  =  8px .  □ 
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Observation:  If  we  choose  /?  to  be  a  suitable  constant 
between  between  0  and  1,  the  above  proof  gives  63  = 
0(y/i 2)  =  0(Vi  +  <5)  though  the  radius  of  the  empty 
ball  becomes  a  constant  fraction  of  px .  Also,  the  entire 
proof  remains  valid  when  we  replace  px  with  f(x). 

4.  Normal  approximation 

We  will  approximate  the  normals  by  the  vectors  from  the 
sample  points  toward  the  centers  of  the  Delaunay  balls 
incident  to  them.  First,  we  derive  an  upper  bound  on 
this  normal  approximation  error  in  Theorem  1.  Then, 
we  describe  a  simple  algorithm  for  approximating  the 
normals  whose  justification  is  given  by  the  theorem  and 
Lemma  4. 

4.1  Analysis 

The  idea  is  to  measure  how  much  one  can  tilt  an  empty 
ball  anchored  at  a  point  p  with  its  center  in  the  direction 
of  Up  while  keeping  it  empty.  The  amount  of  tilt  depends 
on  how  close  the  sample  points  are  and  also  how  big  the 
ball  is. 

Theorem  1.  Letp  G  P  be  incident  to  an  empty  ball  Bc,r 
where  r  >  A  f(p)  and  c  G  flj.  Then, 

sin  (Zpc, up)  <  2(y  +  l)(ei  +  2 Vd)  +  0(5)  +  0(e2) 
A 

for  a  sufficiently  small  s  >  0  and  5  >  0. 

Proof  Let  B  =  Bc,r.  Assume  that  pc  makes  an  angle 
P  with  the  normal  n^.  Let  Bin  and  Bout  be  two  balls 
with  radius  f(p)  and  tangentially  meeting  the  surface  at 
point  p  as  Figure  5  shows.  Let  m  be  the  center  of  Bout. 
We  know  the  surface  £  is  outside  these  two  balls.  By 
Lemma  3  the  ball  B1  =  B(m ,  (1  —  3 5)f(p)),  a  shrunk 
ball  of  Bout ,  is  also  empty  of  samples.  Therefore  no 
sample  point  is  inside  the  shaded  area,  see  Figure  5. 


Figure  5 :  Illustration  for  Theorem  1 . 

Observation  A.  As  p  increases,  the  radius  of  the  inter¬ 
section  circle  of  the  boundaries  of  B  and  B '  increases, 


i.e.  \\w  —  q\\  increases,  and  vice  versa. 

Observation  B.  Suppose  that  when  \\w  —  q\\  = 

s/2sif(p)  the  angle  j3  has  the  claimed  bound.  Then,  if 
we  show  \\w  —  (/||  <  s/2sif(p),  we  are  done  following 
Observation  A. 

Assume  \\w  —  q\\  =  s/2sif(p).  Let  z  be  the  intersec¬ 
tion  point  between  £  and  the  segment  me. 

Consider  the  triangle  formed  by  p,  m  and  c.  We  have 

(1  -  5)f(p)  <  \\m  -  p\\  <  (1  +  5)f(p), 

||c-p||  =  ||c-  w||  >  A f{p) 

and  also  ||c  —  m\\  equals 

s/\\c  —  w\\ 2  —  ||  w  —  q\\ 2  +  s/\\ m  —  w\\ 2  —  \\w  —  q\\ 2. 

If  P  =  Zpc,  up, 

„  lie  —  mil2  —  lie  —  pll2  —  lira  —  pll2 

cos  P  —  - —77 — - — nr - ri - — 

2\\c  -  p\\\\m  -  p\\ 

^  ( VA2  -  24  +  V(1  -  3^)2  -  2£?)2  -  A2  -  (1  +  <5)2 
2A(1  +  <5) 

>  1  -  2(1  +  L2(£i  +  2 VS)2  +  0{52)  +  0{e\). 

A 

Hence 

sin  p  <  VW  —  cos  p) 

<  2(1  +  —)(si  +  2s/5)  +  0(5)  +  0(e\). 

A 

Now  we  show  that  \\w  —  q\\  <  s/2sif(p)  as  required 
by  Observation  B. 

Again,  first  assume  that  \\w  —  q\\  =  y/2eif(p).  We 
can  show  ||p  —  z\\  <  3||p  —  m||  tan p.  Therefore,  from 
the  bound  on  p  ||p  —  ^||  =  0(s\  +  2 V5)f(p)  which  by 
Lipschitz  property  gives  f(z)  <  \/2 f(p)  given  a  suf¬ 
ficiently  small  5  and  e.  We  know  Bz  =  B(z,S\f(z)) 
with  radius  S\f(z)  <  s/2s\f(p)  has  to  contain  at  least 
one  sample  point  by  Lemma  1 .  This  is  impossible  since 
Bz  has  a  radius  at  most  s/2s\f(p)  =  \\w  —  q\\  which 
means  it  lies  completely  in  the  shaded  area.  Therefore, 
llw  —  g||  7^  s/2sif(p).  Now  consider  increasing  Hw  —  ^|| 
beyond  this  distance  while  keeping  z  fixed.  Notice  that 
now  z  is  not  the  intersection  point  between  £  and  the 
segment  me.  It  is  obvious  that  Bz  remains  inside  the 
shaded  area.  Therefore,  again  we  reach  a  contradiction 
to  Lemma  1.  Hence  \\w  —  q\\  can  not  be  larger  than 
V%£l  f(p).  □ 

Implications:  Theorem  1  gives  a  general  form  of 
the  normal  approximation  under  a  fairly  general  sam¬ 
pling  assumption.  One  can  derive  different  normal  ap¬ 
proximation  bounds  under  different  sampling  assump¬ 
tions  from  this  general  result.  For  example,  if  P  is  a 
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(s,  e2,  — ) -sample  we  get  an  0(e)  bound  on  the  normal 
approximation  error.  In  case  P  is  a  (e,  6,  —) -sample,  we 
get  an  0{y/e)  error  bound.  Another  important  implica¬ 
tion  is  that  Delaunay  balls  need  not  be  too  big  to  give 
good  normal  estimates.  One  can  observe  that  if  A  is  only 
y^ma x{£,  5},  we  get  0(52  +  £4)  error.  Algorithmic 
implication  of  this  fact  is  that  a  lot  of  sample  points  can 
qualify  for  normal  estimation. 

Observe  that  the  proof  of  Theorem  1  remains  valid 
even  if  the  sample  point  p  is  replaced  with  any  point 
x  G  M3  meeting  the  conditions  as  stated  in  the  corol¬ 
lary  below.  We  use  this  fact  later  in  feature  estimation. 

Corollary  5.  Let  x  G  M3  be  any  point  with  \\x  —  x\\  < 
dpx  and  Bc^r  be  any  empty  ball  incident  to  x  so  that 
r  =  Ll(px).  Then,  Zxc,  =  0(e  +  s/5)  for  sufficiently 
small  s  and  6. 

4.2  Algorithm 

We  know  from  Theorem  1  that  if  there  is  a  big  Delaunay 
ball  incident  to  a  sample  point  p ,  then  the  vector  from 
p  to  the  center  of  the  ball  estimates  the  normal  direction 
at  the  point  p.  On  the  other  hand,  the  observation  after 
the  proof  of  Lemma  4  assures  that  for  each  point  x  G  E, 
there  is  a  sample  point p  within  0{y/e  +  S)f(x)  distance 
with  an  empty  ball  of  radius  0  (/(#)).  This  means  there 
is  a  big  Delaunay  ball  incident  to  p  where  the  vector  pc 
approximates  and  hence  .  Algorithmically  we  can 
exploit  this  fact  by  picking  up  sample  points  that  are  in¬ 
cident  to  big  Delaunay  balls  only  if  we  have  a  scale  to 
measure  ‘big’  Delaunay  balls.  For  this  we  assume  the 
third  condition  in  the  sampling  which  says  that  the  sam¬ 
ple  is  locally  uniform. 

Let  A p  be  the  distance  of  p  to  its  Kth  nearest  neigh¬ 
bor.  By  sampling  condition  Xp  >  sf(p).  Therefore,  any 
Delaunay  ball  incident  to  p  with  radius  more  than  rXp 
will  give  a  normal  estimation  with  an  error  6(^  +  1) (e) 
according  to  Theorem  1  under  the  assumption  that  P  is 
a  (£,£2,  k) -sample.  It  is  important  that  Xp  is  not  arbi¬ 
trarily  large  since  then  no  Delaunay  may  qualify  for  the 
size  threshold.  This  concern  is  alleviated  by  the  fact  that 
<  e'fip)  where  s'  =  (s  +  6  t8]- 

Notice  that  the  error  decreases  as  r  increases.  How¬ 
ever,  as  we  indicated  before  refip ),  the  radius  of  the 
big  Delaunay  ball,  can  be  as  small  as  52  f(p)  to  give  an 
0(y/e)  error.  This  explains  why  a  large  number  of  De¬ 
launay  balls  give  good  normal  estimations  as  Figure  1 
illustrates. 

ApproximateNormal(P,  t) 

Compute  Del  P\ 

for  each  p  G  P  compute  Xp ; 

i  f  there  is  a  Delaunay  ball  incident  to  p 


with  radius  larger  than  rXp 

Compute  the  largest  Delaunay  ball  Bc,r 

incident  to  p\ 

Approximate  the  normal  direction  at  p  by  pc. 

endif 

Notice  that,  alternatively  we  could  have  eliminated  the 
parameter  r  in  the  algorithm  by  looking  for  the  largest 
Delaunay  ball  incident  to  a  set  of  fc-nearest  neighbors  of 
p  for  some  suitable  k.  Again,  thanks  to  Lemma  4,  we  are 
assured  that  for  a  suitable  k ,  one  or  more  neighbors  have 
Delaunay  balls  with  radius  almost  equal  to  the  medial 
balls.  However,  this  approach  limits  the  number  of  sam¬ 
ple  points  where  the  normals  are  estimated.  Because  of 
our  earlier  observation,  the  normals  can  be  estimated  at 
more  points  where  the  Delaunay  ball  is  big  but  not  nec¬ 
essarily  as  big  as  the  medial  balls.  In  contrast,  as  we  see 
next,  feature  estimation  needs  the  Delaunay  balls  almost 
as  big  as  the  medial  ones. 

5.  Feature  approximation 

We  approximate  the  local  feature  size  at  a  sample  point 
p  by  first  approximating  the  medial  axis  with  a  set  of 
discrete  points  and  then  measuring  the  distance  of  p  from 
this  set.  We  are  guaranteed  by  Lemma  4  that  there  are 
many  sample  points  which  are  incident  to  big  Delaunay 
balls.  The  furthest  Voronoi  vertices  from  these  sample 
points  in  Tlj  and  TIq  approximate  the  inner  and  outer 
medial  axis  respectively.  For  a  point  p  E  P,  we  call  the 
furthest  Voronoi  vertex  from  p  in  Vp  D  Qj  as  the  inner 
pole  p+  of  p.  Similarly  one  may  define  the  outer  pole 
p~  of  p  which  resides  in  Qq. 

In  line  with  the  previous  results  on  medial  axis  ap¬ 
proximation  [4,  6,  7],  we  claim  that  a  certain  subset  of 
the  medial  axis  is  approximated  by  poles.  Let  x  and  x ' 
be  two  points  where  the  medial  ball  B  centered  at  m 
meets  S.  Call  Zxmx 1  the  medial  angle  at  m  if  it  is  the 
largest  angle  less  than  7r  made  by  any  two  such  points  of 
B  flS.  Let  Ma  C  M  be  the  subset  where  each  point 
m  G  Ma  has  a  medial  angle  at  least  a. 

5.1  Analysis 

We  show  that  each  medial  axis  point  mx  with  a  large 
enough  medial  angle  is  approximated  by  a  pole.  The 
idea  is  as  follows.  Consider  the  large  ball  incident  to  a 
sample  point  p  guaranteed  by  Lemma  4.  Then  we  de¬ 
form  it  to  a  large  Delaunay  ball  with  the  center  at  p+ . 
First,  during  this  deformation  the  ball  cannot  be  tilted 
too  much  since  the  vector  from  p  to  the  center  has  to 
approximate  the  normal  by  Theorem  1 .  Second,  the 
center  in  the  tilted  direction  cannot  move  too  much  due 
to  Lemma  6  as  stated  below.  The  result  of  these  con¬ 
straints  is  that  the  center  p+  of  the  Delaunay  ball  remains 
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close  to  the  center  of  the  original  ball  which  in  turn  is 
close  to  mx . 

Lemma  6.  Let  B  =  Bc^r  be  an  empty  ball  whose  bound¬ 
ary  passes  through  a  sample  point  p.  Let  z  be  a  point  on 
E  whose  distance  to  the  boundary  of  B  is  less  than  s' pz 
for  s'  <  1.  Let  B 1  —  Bc>y  be  an  empty  ball  obtained 
by  expanding  B  while  keeping  d  on  the  ray  pc  and  p  on 
its  boundary.  If  r  >  j3pz,  then  we  have 

||c_c,||<  (£i  +  £,)(2/3  +  d  —  £i) 

—  2/3(1  -  cos  Zpcz)  -  2s i  -  2 s'  cos  Zpcz^z' 

Proof  Let  y  be  the  closest  point  to  z  on  the  boundary  of 
B.  Obviously  y ,  c  and  z  are  collinear.  See  Figure  6.  We 
have  \[y  —  z\\  <  s' pz.  Let  z'  be  the  point  where  the  line 
of  c!z  intersects  the  boundary  of  B1 .  Since  a  ball  cen¬ 
tered  at  z  with  radius  Sif(z)  can  not  be  empty  of  sample 
points  by  Lemma  1,  we  have  || z'  —  z ||  <  Sif(z)  <  s± pz. 

Consider  the  triangle  made  by  c,  d  and  z.  For  conve¬ 
nience  set  Lpcz  —  a,  \\c  —  d\\  =  Ac,  \\z  —  z*\\  =  A z 
and  \\y  —  z\\  =  Ay.  We  know 

|| d  —  z\\  =  (Ac)2  +  ||c  —  z\\2  +  2Ac||c  —  z\\  cosa. 

This  gives 

(r+Ac—Az)2  =  (Ac)2  +  (r+A?/)2+2Ac(r+A^)  cosa 

from  which  we  get 

Ac  =  (r  +  Ay)2  -  (r  -  Az)2 

2 (r  -  Az)  -  2 (r  +  Ay)  cosa 
<  (£i  +  d)(2(3  +  s'  —  Si) 

~  2/3(1  —  cosa)  —  2s  i  —  2  s1  cos  a^z’ 

by  plugging  in  A^  <  s1pz,  Ay  <  s'pz  and  r  >  (3pz. 

□ 


Figure  6:  Illustration  for  Lemma  6. 


Theorem  2.  For  each  point  mx  G  Ma  f|  fli  where  a  = 
£4  +5  4,  there  is  a  sample  point  p  within  0(s  *  +  5  *)px 
distance  of  x  so  that  the  pole  lies  within  0(s  *  + 
5*)px  distance  from  mx  where  s  and  5  are  sufficiently 
small. 

Proof.  Consider  the  ball  B  =  Bc^r  guaranteed  by 
Lemma  4  whose  boundary  passes  through  a  sample 
point  p.  Its  radius  r  is  (1  —  2^Jsf)px,  \\p  —  x\\  <  s%px 
and  ||c  —  mx\\  <  2<y/s^px.  Let  B '  =  Bp+^r>  where p+ 
is  the  inner  pole  of  p  and  d  —  \\p  —  p+||.  The  ball  B1  is 
Delaunay  and  its  radius  d  >  r  >  (1  —  2 yjsf)px. 

Focus  on  the  two  balls  B  and  B1  passing  through  p , 
see  Figure  7.  The  ball  B  has  mx  inside  it  which  means 
that  its  radius  is  at  least  (1  —  5)f(p)/2.  So,  the  ra¬ 
dius  of  B '  being  bigger  than  that  of  B  is  also  at  least 
(1  —  5)f(p)/2.  Therefore,  by  Theorem  1  the  vectors  pc 
and  pp+  make  0(e  +  y/S)  angle  with  and  at  most 
double  of  this  angle  among  them.  Let  d  be  the  point  on 
the  segment  pp+  so  that  pc'  has  the  same  length  as  pc. 
Clearly, 

|| C  -  c'||  <  ||p  -  c||Zcpc'  <  (1  -  2 y/e^)0(e  +  VS)px. 

Now  we  can  bound  the  distance  ||c  —  p+||  if  we  have  a 
bound  on  ||c'  —  p+ 1|.  We  will  apply  Lemma  6  to  the  ball 
B"  =  ^cMIp-c'H  to  bound  ||c'  —  p+||.  Since  mx  G  Ma 
there  are  two  points  x  and  x1  in  E  so  that  Zxmxxf  >  a. 
Take  z  in  Lemma  6  as  the  point  x  or  x1  which  makes  the 
angle  Zzmxp  at  least  a/2. 

With  this  set  up  we  show  that  /3  and  s'  in  Lemma  6 
are  1  —  0(y/s  +  S)  and  0(s  +  y/S)  respectively.  Since 
the  radius  of  Bn  is  r  >  (1  —  2-sfsf)px  —  (1  —  2-sfsf)pz, 
the  claim  for  /3  follows. 


Figure  7:  Illustration  for  Theorem  2.  The  ball  Bc^r  is 
deformed  to  the  Delaunay  ball  B1  =  Bp+  y.  The  ball 
B "  =  Bci^p_c>\\  on  the  right  is  a  shrunk  version  of  B1 . 

For  s’ ,  consider  the  point  y  where  the  ray  dz  meets  the 
boundary  of  Bn ,  refer  to  Figure  7.  We  have  |  |cr  —  mx  \  |  < 
||c  —  mx  ||  +  ||c  —  c7 1|  =  0(y/s  +  y/~5)px  and  hence 
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\\y-z\\  =  \W  -  z\\  -  \\cf  -  y\\ 

<  Urns  -  2: 1 1  +  ||c'  -  ra#J|  -  He'  -  y\\ 

<  Pz  +  0(V^  +  y/$)pz  —  (1  —  ^\f^l)pz 

—  0{y/e  +  V6)pz. 

So,  we  can  apply  Lemma  6  with  s'  =  0(y/e  +  \/£),  and 
/3  =  1  —  0(yjs  +  5).  Observe  that,  since  the  points  c' 
and  ra^  are  nearby,  the  angle  Zpc'y  is  almost  equal  to 
Zzmxp.  So,  we  can  safely  assume  Zpc'y  >  f.  With 
a  =  ei  +  5v  Lemma  6  gives 

||p+  — c'll  =  (0(£^  +  +  di))pz 

=  0{e*+8*)pz. 

The  claim  of  the  theorem  follows  as 

\\p+-mx\\  <  \\p+  -c'||  +  ||c'  -mx\\ 

=  0(e*  +Si)px  +  0(ei  +  5*)px 
=  0(e*+6*)px. 

□ 

For  each  point  x  G  S  where  ra^  G  Ma  the  previ¬ 
ous  theorem  guarantees  the  existence  of  a  sample  point  p 
whose  pole  approximates  mx .  Actually,  the  proof  tech¬ 
nique  can  also  be  used  to  show  that  any  Delaunay  ball 
with  radius  almost  as  big  as  px  and  incident  to  a  sample 
point  close  to  x  has  its  center  close  to  mx . 

Theorem  3.  Let  x  G  S  be  a  point  so  that  mx  G  Ma 
for  a  =  £4  +  4 .  Then  for  any  point  p  G  P  within 

£3 px  distance  of  x  and  with  an  incident  Delaunay  ball 
of  radius  at  least  (1  —  0{\fe  +  s/S)px,  the  pole  p+  lies 
within  0(£8  +  5*)px  distance  from  mx. 

sketch.  Notice  that  if  and  make  small  angle,  we 
will  be  done.  Then,  we  have  two  segments  pp+  and  xmx 
almost  parallel  where  p  and  x  are  close.  Also,  these  seg¬ 
ments  can  be  shown  to  be  of  almost  same  lengths  by  the 
given  condition  and  a  proof  similar  to  that  of  Theorem  2. 
This  would  imply  mx  and  p+  are  close. 

Observe  that  we  cannot  assert  that  is  small 

directly  from  the  normal  variation  lemma  in  Amenta  and 
Bern  [2]  as  the  distance  among  them  is  at  most  £3 px  by 
Lemma  4  (notice  px  vs.  f(x)).  Since  p  and  x  are  at 
most  £3 px  apart  and  the  distance  of  to  p  and  hence 
to  x  is  D(px),  App+x  =  0(£ 3).  By  Corollary  5,  it 
can  be  shown  that  p+x  makes  0(y/i 3)  angle  with  n^. 
Therefore,  pp+  makes  3)  angle  with  n^.  It  is  easy 
to  show  that  px  is  at  least  Ll(f(p)).  So,  the  angle  be¬ 
tween  pp+  and  Up  is  0(s/£3)  =  0(6  8  +5 s).  One 
can  deduce  that  1 1 p+  —  mx  \  |  has  the  claimed  bound  with 
Znp,nx  =  0(^).  □ 


5.2  Algorithm 

Theorem  2  and  Theorem  3  suggest  the  following  algo¬ 
rithm  for  feature  estimation  at  any  point  x  G  S  where 
mx  G  M  1  1 .  Theorem  2  says  that  x  has  a  sample 

£4-|-<)4 

point  p  within  a  neighborhood  of  £3  px  whose  pole  p+ 
approximates  mx .  Also,  Theorem  3  says  that  all  sam¬ 
ple  points  within  £3  px  neighborhood  of  x  with  a  large 
enough  Delaunay  ball  have  their  poles  approximate  mx. 
Therefore,  if  we  take  the  pole  of  a  sample  point  q  whose 
distance  to  q  is  largest  among  all  sample  points  within  a 
neighborhood  of  x,  we  will  get  an  approximation  of  mx . 

We  search  the  neighborhood  of  x  by  taking  k  nearest 
neighbors  of  a  sample  point  s  close  to  x.  If  we  assume 
that  P  is  a  (5, 5:  k) -sample  for  some  k  >  0,  ^-nearest 
neighbors  cannot  be  arbitrarily  close  to  x.  Notice  that  if 
we  do  not  prevent  oversampling  by  the  third  condition 
of  noisy  sampling,  we  cannot  make  this  assertion.  In  the 
algorithm,  we  simply  allow  an  user  supplied  parameter 
k  to  search  the  k  nearest  neighbors.  Since  we  want  to 
cover  all  points  of  E,  we  simply  take  all  points  of  P  and 
carry  out  the  following  computations. 

For  each  point  p  G  P  we  select  ^-nearest  neighbors 
for  a  suitable  k.  Let  Np  be  this  set  of  neighbors.  First, 
for  each  q  G  Np,  we  determine  the  Voronoi  vertex  vq 
in  Vq  which  is  furthest  from  q.  This  is  one  of  the  poles 
of  q.  Let  (p)  =  || vq  —  q\\.  Select  the  point  pi  G  Np 
so  that  l\{p\)  is  maximum  among  all  points  in  Np.  By 
Theorem  2  and  Theorem  3,  vPl  approximates  a  medial 
axis  point  mx  if  x  G  M  1  1 .  However,  we  do  not 

£4  4 

know  if  mx  is  an  inner  medial  axis  point  or  an  outer  one. 
Without  loss  of  generality  assume  that  mx  is  an  inner 
medial  axis  point.  To  approximate  the  outer  medial  axis 
point  for  x ,  we  determine  the  Voronoi  vertex  uq  in  Vq  for 
each  q  e  Np  so  that  quq  makes  more  than  |  angle  with 
PiVPl.  Let  I2 (q)  =  || uq  —  q\\.  Then,  we  select  the  point 
P2  £  Np  so  that  €2(^2)  is  maximum  among  all  points  in 
Np.  Again,  appealing  to  Theorem  2  and  Theorem  3  for 
outer  medial  axis,  we  can  assert  that  uP2  approximates  a 
medial  axis  point  for  x. 

ApproximateFeature(P,  k) 

Compute  Del  P\L  f  \ 

for  each  p  G  P  compute  k  nearest  neighbors  Np\ 
compute  pi  G  Np  whose  distance  to 
its  pole  vPl  is  maximum 
among  all  points  in  Np\ 
compute  P2  £  Np  with  a  pole  v P2  so  that 
^P2Vp2,PivPl  >  f  and  \\p2  - 
is  maximum  among  Np\ 

L  :=  LU  {vPl,Vjto}\ 
endf or 

for  each  p  G  P  compute  the  distance  of  p  to  L. 

As  we  have  observed  already,  a  subset  of  the  medial 
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Figure  8:  Left:  Medial  axis  approximated  by  centers  of  big  Delaunay  balls  for  a  noisy  HORSE.  For  a  chosen  threshold, 
some  parts  of  the  legs  do  not  have  medial  axis  approximated  though  still  many  centers  lie  near  the  surface.  Right: 
Medial  axis  well  approximated  by  the  poles  as  computed  by  ApproximateFeature. 


axis  is  not  approximated  by  the  poles.  These  are  exactly 
the  points  on  the  medial  axis  which  have  a  small  medial 
angle.  This  type  of  exclusions  are  also  present  in  earlier 
medial  axis  approximation  results  [4,  6,  7].  The  implica¬ 
tion  of  this  exclusion  is  that  features  cannot  be  properly 
estimated  for  points  whose  closest  point  on  the  medial 
axis  resides  in  the  excluded  part.  However,  if  the  sam¬ 
pling  is  sufficiently  dense,  the  excluded  part  is  indeed 
small  in  most  cases.  Figure  8  shows  the  result  of  feature 
approximations  for  a  model  in  three  dimensions. 


References 

[1]  M.  Alexa,  J.  Behr,  D.  Cohen-Or,  S.  Fleishman,  D. 
Levin  and  C.  Silva.  Point  set  surfaces.  Proc.  IEEE 
Visualization  (2001),  21-28. 

[2]  N.  Amenta  and  M.  Bern.  Surface  reconstruction  by 
Voronoi  filtering.  Discr.  Comput.  Geom.  22  (1999), 
481-504. 

[3]  N.  Amenta,  M.  Bern  and  D.  Eppstein.  The  crust 
and  the  /?- skeleton:  combinatorial  curve  recon¬ 
struction.  Graphical  Models  and  Image  Process¬ 
ing ,  60  (1998),  125-135. 

[4]  N.  Amenta,  S.  Choi  and  R.  K.  Kolluri.  The  power 
crust,  union  of  balls,  and  the  medial  axis  transform. 
Comput.  Geom.:  Theory  Applications  19  (2001), 
127-153. 

[5]  D.  Attali  and  A.  Montanvert.  Modeling  noise  for 
a  better  simplification  of  skeletons.  Proc.  Internal 
Conf.  Image  Process.  3,  13-16,  1996. 


[6]  J.  D.  Boissonnat  and  F.  Cazals.  Natural  neighbor 
coordinates  of  points  on  a  surface.  Comput.  Geom. 
Theory  Applications  (2001),  87-120. 

[7]  F.  Chazal  and  A.  Lieutier.  Stability  and  homotopy 
of  a  subset  of  the  medial  axis.  Proc.  ACM  Sympos. 
Solid  Modeling  and  Applications  (2004),  243-248. 

[8]  T.  K.  Dey  and  S.  Goswami.  Provable  surface  re¬ 
construction  from  noisy  samples.  Proc.  20th  Annu. 
Sympos.  Comput.  Geom.  (2004),  330-339. 

[9]  T.  K.  Dey  and  J.  Sun.  Adaptive  MLS  surfaces  for 
reconstruction  with  guarantees.  Proc.  Eurograph¬ 
ics  Sympos.  Geom.  Processing  (2005),  43-52. 

[10]  T.  K.  Dey  and  W.  Zhao.  Approximating  the  medial 
axis  from  the  Voronoi  diagram  with  a  convergence 
guarantee.  Algorithmica  38  (2003),  179-200. 

[11]  R.  Kolluri.  Provably  good  moving  least  squares. 
Proc.  ACM- SIAM  Sympos.  Discrete  Algorithms 
(2005),  1008-1017. 

[12]  B.  Mederos,  N.  Amenta,  L.  Velho  and  H.  de 
Figueiredo.  Surface  reconstruction  from  noisy 
point  clouds.  Proc.  Eurographics  Sympos.  Geom. 
Processing  (2005),  53-62. 

[13]  N.  Mitra,  J.  Nguyen  and  L.  Guibas.  Estimating  sur¬ 
face  normals  in  noisy  point  cloud  data.  Internat.  J. 
Comput.  Geom.  &  Applications  (2004). 

[14]  M.  Pauly,  R.  Keiser,  L.  Kobbelt  and  M.  Gross. 
Shape  modeling  with  point- sampled  geometry. 
Proc.  ACM  SIGGRAPH  2003 ,  641-650. 

[15]  www.  c se .  ohio-  state .  edu/~ tamaldey /normfet . html 


10 


